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Abstract 

Using a known result by Falk et al. for the 1/ml correction to the dilepton invariant 
mass spectrum in the decay B — ► X s ^ + fi^, we calculate the l/rn^ correction to the 
left-right muon polarization asymmetry in this decay. Employing an up-to-date range of 
values for the non-perturbative parameter Ai, we find that the correction is much smaller 
than it should have been expected from the previous work by Falk et al. 
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Introduction 



Rare decays of B mesons have been studied extensively in the last few years. Such decays 
are forbidden in the tree-level approximation and proceed through loop diagramms only. 
Consequently, they are sensitive to the complete particle content of a given theory and their 
analysis may thus shed some light on possible new physics beyond the standard model (SM). 
In this paper, we focus on the inclusive decay B — ► X s fi + n~ , where X s denotes an arbitrary 
hadronic final state with total strangeness —10- in contrast to the decay B — > X s 7, the three- 
body decay B — > X s allows one to define and measure several kinematic distributions. 
Beside the invariant mass spectrum of the lepton pair a forward-backward charge asymmetry 

and a left-right polarization asymmetry [Q, || have been proposed. The simultaneous 
measurement of these distributions allows one to extract the values of the Wilson coeffients 
which govern the decay amplitude and contain the short distance physics M. The values of 
these Wilson coefficients are sensitive to new physics and in a measurement of the kinematic 
distributions one might be able to detect deviations from the SM prediction. For this reason 
the distributions mentioned above should be calculated as precise as possible. 

Meanwhile a complete next-to-leading order (NLO) calculation of the relevant Wilson 
coefficients is available |5|,||. These coefficients properly include the short distance QCD 
effects and the NLO approximation considerably reduces the (otherwise significant) theoretical 
uncertainty in the final result which stems from its dependence on the renormalization scale. 

Apart from these perturbative calculations which rely entirely on the spectator model 
(defined as the free quark decay model including perturbative QCD), long distance effects 
due to the non-perturbative physics at low scales have to be taken into account. Important 
corrections are first of all generated by intermediate cc bound states (e. g. the channel b — > 
s J/ip — > which influence the kinematic distributions even far away from the resonance 

region. For a way to deal with these effects we refer to the literature [7-1C]. 



Additionally, there are non-perturbative effects due to the binding of the b quark inside 
the B meson. It is strongly expected that for inclusive decays of heavy quarks such effects 
may be incorporated by a controlled expansion in inverse powers of the mass of the decaying 
quark, the so-called heavy quark expansion (HQE) pd| , [T^| p^|. The first term of this expansion 



can be shown to reproduce the spectator model. Falk, Luke and Savage [14| calculated the 
first non-vanishing - i.e. 0(l/m|) - correction to the dilepton invariant mass spectrum in 
the decay B — > X s some time ago. 

It is this type of corrections we are dealing with in the following. We slightly extend the 
work of Liu and Delbourgo |J about the left-right polarization asymmetry by presenting the 
1/m^ correction to this distribution. It should particularly be stressed that in the numercial 
analysis we take into account a present range of values for Ai (one of the non-perturbative 
parameters that occur by performing the HQE), and that these values differ drastically from 
the one used in Ref. f\M. 



Since we neglect the mass of the leptons, all statements in this paper apply likewise to the decay B — » 
X a e + e~. However, this decay mode is experimentally even harder accessible than B — * X s [i + fj,~, so we 
decided to mention only the latter one explicitly in the text. 

2 - ) Whether this expansion can be justified theoretically rigorous is still the subject of ongoing discussions [ p^[ . 
At the moment, phenomenology seems to be the only way to judge the validity of the HQE. 
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The left-right asymmetry 

The general framework for decays like B — * X s /i + ' \x~ is the effective theory approach with the 
effective Hamiltonian as the central element. Once the effective Hamiltonian is given one can 
calculate the decay amplitude and the kinematic distributions mentioned above. The effective 
Hamiltonian for the decay B — > X s has been calculated in the NLO approximation by 

Misiak || and independently by Buras and Miinz ||. In our paper we entirely use the 
notation introduced in Ref. Q and we refer to this paper for the operator basis, the analytic 
expressions of the Wilson coefficients and for the details concerning the NLO calculation. 
Let us now consider the left-right asymmetry. It is defined as 



dA LR 
ds 



dT L dT 



R 



ds 



ds 



where 



mt 



is the scaled 4- momentum transfer to the muons and dT L /ds (dT R /ds) denotes the invariant 
mass spectrum for a decay into purely left-handed (right-handed) muons. Instead of calculat- 
ing dT L /ds and dT R /ds separately, dA LR /ds can alternatively be obtained directly from the 
invariant mass spectrum dT /ds. With the operator basis used in Ref. ||, dT/ds is derived 
in terms of the Wilson coefficients Cg and C\q. The corresponding operator Qg contains a 
vector coupling of the muons, Q±o an axial-vector coupling. Defining a left-handed and a 
right-handed operator by 



Ql 
Qr 



Q9 — Q10, 
Q9 + QiO) 



one obtains dT/ds in terms of the corresponding Wilson coefficients Cl and Cr by the simple 
substitution 

Cg = Cl + Cr, 
C\o = Cr - Cl- 

Since the remaining operators Q\ 
in this basis 



Qs contain only vector couplings to the muons, one finds 



dA LR 
ds 



dT 

ds 



c R =o 



dT 

ds 



C L =0 



The final result may be transformed back to the basis Qg, Qiq and re-expressed in terms of 
Cg and Cio0- Following these steps one reproduces the result for the left-right asymmetry 
in Ref. || if one starts with the invariant mass spectrum as given in Ref. 0. Keeping this 



3) Note that 



dT 

ds 



dT L 
ds 



and — 

ds 



dr R 

ds 



because t5 contains right-handed components due to the operators Qi 

" s I C n~ 

same terms present in | c _ Q by taking the difference. 



■ Qs which cancel against the 
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procedure in mind one can extract the l/m 2 -correction to the left-right asymmetry from the 
1/mjj-correction to the invariant mass spectrum given in [14|. We stress that it is not necessary 



to repeat the full calculation Falk et al. have done in order to find the l/m 2 -correction for 
the left-right asymmetry^]. We obtain 



dA LR 

ds 
dA^ R 

ds 
dAr/ R 

ds 



dA^ R 
ds 
T(b 
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dA\ R 

ds 
cl~V) 



m 

T(b- 



\V t 



K{Z 



cl v) 



tb\ 



Vf 



where we defined z 



fcf*(Ai,l) 



/(*) 
= m c /mb and 



tb\ 



v cb 

Vts 



Vcb 



a 

■2 



1 - s) 2 C w [(1 + 2s)Re ((7 9 e// ) + 6 Of 
£^ (1 " s)Cio (k% R (Xi,s) + k% R (\ 2 , s 



(0)e// 



2ml 



i_ { (-2§ 2 + 3s + 5) Re (c e g ff ) - 2(3S - 5)cf e// } , 
{3 (-10s 2 + 15s + l) Re ((7 9 e// ) - 6(7s - 5)cJ 0)e// } . 



dAQ R /ds denotes the spectator model result. It agrees with the result in Ref. |J. The func- 
tions f(z) and k{z) are the phase-space factor and the one gluon correction to the decay b — > 



ceP given in Eqn. (2.31) and (2.32) of Ref . §. The Wilson coefficients Cf ;e// , Cg^and Ci 
are also defined in Ref. g, Eqn. (2.3), (2.28) and (2.8). Note that the effective coefficient 
Cg^ depends explicitely on s because it contains the influence of the operators Qi • ■ Qe- 
Finally, dA^/ds denotes the 1/m^-correction, given in terms of the two HQET parameter 



Ai and A2. A definition of them as matrix elements can be found in Ref. [14], Eq. (2.21) 



Numerics 

Let us now examine the size of the correction compared with the spectator model result. In 
doing this we do not incorporate the long distance effects due to intermediate cc resonances. 
We do not neglect them because they are small - in fact, as already mentioned, they do 
have a large effect on the shape of the spectra even far away from the resonance region. The 
common way to deal with these resonances is a simple replacement of by Cg^ + Res(s) 
(see e.g. Ref. |^0|), where the function Res(s) is chosen to produce resonance peaks at the 
masses of the intermediate cc bound states. The relative correction to the spectator model, 
however, is nearly the same with or without this replacement of Cg . So, for simplicity, we 
show the non-resonant results only. 

In the appendix we give a list of the numerical values of all input parameter we used. 
Some remarks should be made concerning the non-perturbative parameter Ai and A2- From 
the B — B* mass splitting the numerical value of A2 is well known to be 0.12 GeV 2 [p!q[ . Ai is 
much less understood. In Ref. (13] Ai = 0.5 GeV 2 was used in the numerical analysis. These 
values led to an enhancement of roughly 10 percent over the full range of s in the invariant 
mass spectrum. This value of Ai, however, is no longer appropriate. Even if there is still 
some controversary about the magnitude, the sign of Ai is strongly believed to be negative. 
There have been attempts to calculate Ai using QCD sum rules [fl6| — P~9|| and, alternatively, to 

4 'This is not the case for the l/m^-correction to the forward-backward asymmetry, which cannot be extracted 
from the result in [h 411 . 
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extract it by a fit to the data from B and D decays p0|,plj|. Even an upper bound has been 
proposed |2^-p4|. However, the values found lie in the range -0.6 GeV 2 < Ai < -0.1 GeV 2 , 
which is quite large. 

In Fig. 1 we plot the spectator model prediction with and without the 1/m 2 correction. 
The influence of the correction never exceeds a rate of 3.5% in the high s-region (s > 0.4) 



10 



< 4 



2 - 




Figure 1: The left-right asymmetry in the decay B — > X s /x + /x~ (normalized to the total decay width 
of the B meson). The solid line represents the spectator model result. The shaded area corresponds to 
the left-right asymmetry including the 1 /m\ correction, varying Ai between -0.6 GeV 2 and -0.1 GeV 2 . 
The upper edge of the area is given by Ai = —0.1 GeV 2 , the lower edge by Ai = —0.6 GeV 2 . 



and maximally reduces the spectrum by 9.5% for very small values of s. We emphasize that 
the smallness of the correction is due to an accidental cancellation of the corrections k^ R and 
k^- With Ai ~ — 3A2 the functions k^ R and k% R almost cancel in the high s region. For the 
integrated left-right asymmetry the correction lies between +3% and -5% compared to the 
spectator model result. 

This small correction in mind it is interesting to show again the l/m 2 correction to the 



invariant mass distribution with negative values of Ai. This is done in Fig. 2. 
correction is significantly reduced compared to the result in Ref. 



14 



The overall 
As in the case for the 

left-right asymmetry, the correction to the spectator model is at most 3.5% for s > 0.4 and 
maximally differs by 5.2% near s = 0.1. The correction results in a maximal decrease of the 
integrated spectrum of 3.8% for Ai = —0.6 GeV 2 . In comparison, the integrated spectrum 
increases by 9.3% taking Ai as high as the value used in Ref. [14]. 

As long as Ai is poorly known, the following numerical formulae might be useful. They 
describe the relative size of the 1 /m 2 correction to the spectator model results of the integrated 
left-right asymmetry and the branching ratio: 



AA LR 



Ai 
GeV 2 



+ 0.28 x 15.3' 
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Figure 2: The dilepton invariant mass spectrum in the decay B — ► X s /i + fi~ (normalized to the 
total decay width of the B meson). The solid line represents the spectator model result. The shaded 
area corresponds to the spectrum including the 1/to 2 correction, varying Ai between — 0.6 GeV 2 and 
—0.1 GeV 2 . As in Fig. 1 the upper edge of the area is given by Ai = —0.1 GeV 2 , the lower edge by 
Ai = —0.6 GeV 2 . The dashed line is the result obtained by Falk et al. using Ai = 0.5 GeV 2 . 



AB[B -> X s n~yi + 
B[B -> X sf i-fi+}( 



Ai 



GeV^ 



+ 0.28 x 12.0' 



Note that the 1/m 2 correction depends linearly on Ai, therefore these formulae are valid for 
all values of Ai. They also imply that for Ai = —0.28 GeV 2 (and A2 = 0.12 GeV 2 ) one obtains 
the spectator model results 

T£A% R = 4.74 xlO" 6 , 
B[B -» X s fi-n + ]o = 6.03 x 10" 6 . 

It is instructive to compare the size of the 1/m 2 correction to the left-right asymmetry 
with that of other unkown corrections or uncertainties preventing a precise theoretical pre- 
diction of this quantity. We looked at various sources of uncertainties, a complete list of the 
corresponding errors is given in Table 1. 

The most important error stems from the lack of the next-to-next-to-leading order QCD 
calculation. This truncation of the perturbative series manifests itself in the renormalization 
scale dependence of the spectra and integrated rates. By varying the renormalization scale 
[i between and 2m;, (i.e. in such a manner that lnm^/i remains small) we estimated 

the corresponding uncertainty to be about ±8%. Apart from this we have to cope with our 
insufficient knowledge of the masses of the top and charm quark, which enter through loop 
diagrams and the phase-space factor f(z). Both of these masses give rise to an uncertainty of 
roughly ±6%. The errors due to the CKM matrix elements and the semileptonic branching 
ratio are not large, but also comparable to the 1/m 2 correction. The other inaccurate known 
paramters, including a s (Mz) and the masses of the remaining quarks, affect the left-right 
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Uncertainty due to ... 


A A LH 1 A Lti 

AA n /Aq 


Ad/d 


long-distance effects 


±15% 


±10% 


nib/2 < n < 2mb 


±8% 


±2% 


m t = (167 ± 6) GeV 


±6% 


±7% 


m c = (1.4 ±0.1) GeV 


±6% 


±6% 


B[B -» XJ-V] 


±4.8% 


±4.8% 


iVu/Vebf* = 0.95 ± 0.04 


±4.2% 


±4.2% 


a s (M 2 ) = 0.118 ±0.003 


±1.5% 


±0.2% 


m fe = (4.8 ± 0.3) GeV 


±0.7% 


±1.1% 


m s = (0.1 ±0.1) GeV 


±0.7% 


±0.6% 



Table 1: Various sources of uncertainty for the integrated left-right asymmetry as well as for the total 
branching ratio of the decay B — » X s fi + n~ , in order of their numerical importance. For the first two 
entries, see the explanations in the text. 



asymmetry only by a small amount (i. e. not more than ±1%). Apart from these rather 
trivial dependences on the parameters, there remains one potential systematic error, namely 



the long-distance contributions due to cc-resonances. Modelling the resonances as in Rcf. [10|, 
we estimated them to enlarge the left-right asymmetry by roughly 15% (whereby this number 
depends to some extent on the applied experimental cuts). In the worst case, lacking anything 
better, the associated error should be taken to be in the same order of magnitude than 
the correction itself, albeit we would call this a very conservative and perhaps unnecessary 
pessimistic view. 

The numbers of Table 1 show that the l/m| corrections are numerically not very significant 
compared to all the other uncertainties presently existing. However, two remarks should be 
made. On the one hand, we would have arrived at a much larger correction of about 10% 
if we had used the old value Ai = 0.5 GeV 2 as in Ref. |Q (which was published only three 
years ago). It is interesting to note that at the moment most experimental and theoretical 
developments are in favour of the assumption that the l/m 2 , corrections are generically small. 
On the other hand, the l/m 2 corrections are quantities of principle interest. Assuming the 
validity of the HQE, there existence is a fundamental and predictable property of QCD. So 
they should be calculated indepently of their actual size (which, of course, is by no means 
known a priori). Even if today the overall error is dominated by the uncertainties in the SM 
parameters, one cannot say how this picture will change during the forthcoming years. 

Summary 

In this brief note we presented the nonperturbative 1/m 2 correction to the left-right asymme- 
try for the rare decay B — > X s With a present range of values for the parameter Ai we 
found a correction as small as roughly four percent for the integrated left-right asymmetry. 
We also showed the invariant mass spectrum with these new values for Ai and found the 



correction significantly smaller than in Ref. [14]. 
Note added in proof 

After completion of this work a recalculation [25 of the l/m 2 -correction to the decay B — > 
X s fi + fi~ appeared. The authors of Ref. [25] use the same method as Falk, Luke and Savage 
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||14| but obtain a different result. If Ref. pq] should be confirmed, our numerical results would 
of course change. 

Appendix: Input parameter 

m 6 = 4.8GeV, m c = 1.4 GeV, m s = 0.1, = 0, m t = 167 GeV, My^ = 80.2 GeV, 
li = 5 GeV, \V tb \ = 1, \V ts /Vcb\ = 0.95, sin 2 9 W = 0.23, a = a(M z ) = 1/129, a s {M z ) = 
0.118, £(£ -» XJTV) = 0.105, A 2 = 0.12 GeV 2 , -0.6 GeV 2 < Ai < -0.1 GeV 2 (see 
text). 
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